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Quantum decoherence and the transition to semiclassical behavior during inflation has been exten- 
sively considered in the literature. In this paper, we use a simple model to analyze the same process 
in ekpyrosis. Our result is that the quantum to classical transition would not happen during an 
ekpyrotic phase even for superhorizon modes, and therefore the fluctuations cannot be interpreted 
as classical. This implies the prediction of scale-free power spectrum in ekpyrotic/cyclic universe 
model requires more inspection. 

I. INTRODUCTION 

From cosmological observations we know that the current universe is to a good approximation flat, homogeneous 
and isotropic on large scales [U [2] . It is well known that in standard Big Bang cosmology this requires an enormous 
amount of fine-tuning on the initial conditions. Two mechanisms are provided to be possible explanations. The 
first is inflation [5J 0], a period of accelerated expansion occuring between the big bang and nucleosynthesis. The 
second is ekpyrosis [MS], a period of ultra-slow contraction before Big Bang/Big Crunch to an expanding phase. Both 
mechanisms not only manage to address the standard cosmological puzzles but also have the ability to imprint scale 
invariant inhomogeneities on superhorizon scales via a causal mechanism [3l H0H14j . These inhomogeneities are 
thought to provide the seeds which later become the temperature anisotropies in the Cosmic Microwave Background 
and the Large Scale Structure in the universe. This framework of the cosmological perturbation theory is based on 
the quantum mechanics of scalar fields, where the relevant observable is the amplitude of the field's Fourier modes 
[15] . Although treated as a quantum mechanical variable, this amplitude is expected to be stochastic variables, 
characterized by averages of their products, i.e., power spectrum. This interpretation proves to be very accurate in 
the CMB and Large Scale structure analyses. 

However, in order to make this stochastic interpretation consistent, the density matrix has to be diagonal in the 
amplitude basis. This criterion implies that interference terms in the density matrix are highly suppressed and can be 
neglected [TTJ [T5] . Interference is associated with the coherence of the system, i.e., the coherence in the state between 
different points of configuration space [THl HO] ■ A measure of this is the coherence length which gives the configuration 
distance over which off-diagonal terms are correlated [2T] , 

An isolated system described by the Schrodinger equation cannot lose its coherence; a pure state always remains 
pure. However, if it is coarse grained, it may evolve from a pure to a mixed state. One way to realize coarse graining 
is to let the system interact with an environment |19j . The environment consists of all fields whose evolution we are 
not interested in. The state of the system is obtained by tracing over all possible states of the environment. Now, 
even if the state describing system plus environment is pure, the state of the system alone will in general be mixed. 

In the literature, there are various arguments and calculations suggesting that a form of such environment deco- 
herence can indeed occur for inflationary perturbations [211131] ■ The coherence length decreases exponentially for 
wavelengths greater than Hubble radius. Thus perturbations become classical once their wavelength exceeds the Hub- 
ble radius. All of these results lend support to the usual heuristic derivation of the spectrum of density perturbations 
in inflationary models. In this paper, we use a simple model to study whether decoherence can also occur in the 
ekyprotic phase. We find that the coherence lengths continue increasing even for the modes outside the horizon. 
Therefore, the heuristic argument that the modes become classical when they leave the horizon is invalid in the 
ekyprotic phase and requires more careful inspection. 

II. THE MODEL 

A crucial question is how to model the environment. Any realistic model will be very complicated and hard to 
analyze. However, the basic physics should emerge from the simplest models. Hence, we choose a model [21j which 
can be solved exactly: the system is a real massless scalar field <j>\, and the environment is taken to be a second 
massless real scalar field 4>2 interacting with through their gradients. 

The action of system and environment is 

S = J d A xC = f d^x^-g X - (-8^10^1 - ^ 2 d^ 2 - IcB^fa) (1) 
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where g is the determinant of the background metric which is given by 

ds 2 = a 2 {r } ){-dr ] 2 + dx 2 ) (2) 

and c -C 1 is the coupling constant describing the interaction between two fields. Note that this Lagrangian is 
quadratic in the derivative of the fields and can hence be diagonalized for which the interaction term disappears and 
the whole Lagrangian becomes a free field theory. If there is no other field or interaction in our universe, this argument 
is true. However, we suppose there is a hidden interaction such that we can only obeserve the first field <pi but not 
the environment 2 . In other words, we assume the environment and the observed system do not form the diagonal 
basis. This assumption is reasonable since any observed scalar fields (whose reduced density matrix we want) will 
interact with gravitational perturbations (which is a part of the environment). 
Then, the canonical momenta 7Tj conjugate to the fields 0j, i = 1,2 are 

TTj = ^ = a 2 (fa + c0 2 ) (3) 
ofa v ' 

7T 2 = = a 2 (fa + cfa) (4) 
o4> 2 v 7 

where "•" denotes derivative with respect to rj. This allows us to write the Hamiltonian H as 
H = J ePxfcii = J d3x { 2a 2( 1 1 _ C 2) K + ^ 2 - 2c7r l7 r 2 ) + — [(V0i) 2 + (V0 2 ) 2 + 2c(V0i) • (V0 2 )] J (5) 



To study decoherence, it is more convenient to use the functional Schrodinger picture[32 a . The commutation relation 
[0j(x), 7Tj(y)] = iSij5 3 (n — y) is equivalent to making the replacement 7Tj(x) — > — z g^A^ ■ The wave functional 
^[01! 02] obeys the Schrodinger equation 



i~y = HV (6) 
a?7 



We make a Gaussian ansatz for \1/ to be able to find the vacuum or ground state solution: 



* [01,02] = AAexp 



d d xA(0i(x)0 1 (y) + 2 (x)0 2 (y))v4(x,y,? ? ) + 20i(x)0 2 (y)B(x,y,r;) 



(7) 



Note that we have already used the 0i -H> 2 symmetry of the Lagrangian. Furthermore, because of the x f-> y 
symmetry of the above integration, we have to require 

A(x,y,ry) = A(y,x,rj) (8) 
B(x,y,ry) = S(y,x,ry) (9) 

Plug Eq. ([7]) into Schrodinger equation ([6]), it is not difficult to get 

' ' U(X - y -' /) - f (Pz ^,* ^ [A(x, z, 7 ? )A(y, z, ?/ ) + B(x, z, v )B(y, z, 77) - 2cA(x, z, q)fl(y, z, 7/)] 

-f^W(x-y) (10) 



2 dr\ J 2a 2 (l-c 



' 11 X ' ' V ' ' - f d 3 z ^ [B(x, z, r))B(y, z, 77) + A(x, z, » 7 )A(y ) z, 77) - 2cB(x, z, ?y)A(y, z, 77)] 



2 drj J 2a 2 (l-c 



V^(x-y) (11) 



-2cB(x, z, 7 7 ) J B(y ) z, »/)] + y • 2cV 2 5 3 (x - y) (12) 
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.9 In TV 

i — 

drj 2a 2 ( 



^— - J d 3 z[2A(z,z, V )-2B(z,z, V )} 



(13) 



All the above equations come from the comparison of the coefficients in front of 0j(x)0j(y). It is easy to see that 
Eq. pTI and Eq. (10) are equivalent, which is just the result of the symmetry of 0i and 2 . In order to satisfy 



Eq. ( |10[ )-( 12 ), we have to require B(x, y, 77) = cA(x, y, 77), which gives 

¥[0i,02] = Afexp j-^ J d 3 xd 3 y[0i(x)0 1 (y) + 2 (x)0 2 (y) + 2c0 1 (x)0 2 (y)]^(x,y,?7) 



.dln/V 
di] 

^A(x,y,??) 
drj 



^ J d 3 zA{z,z, V ) 

i /" d 3 zA(x, z.^y.z, 77) +a 2 Vj<J 3 (x-y) 



It is more convenient to solve Eq. ( 16 1 in momentum space. Upon writing 



0i(x) 



r0i(k)e' 



(2tt) 3 



we get 



. dA(k, v ) 
drj 



2;„2 



r A 2 (k,7i) -a 2 k 



(14) 

(15) 
(16) 

(17) 
(18) 

(19) 



Here we have already used the relation A(— k, 77) = A(k, 77) coming from Eq. Note that A(k, 77) is only a function 
of |k|, so we will write it as Ak(rj) from now on. This differential equation can be easily solved by assuming 



Uk(v) _ a(rf) 



Then Eq. ( 19 1 becomes 



Uk 



k z --)u k = 
a , 



(20) 



(21) 



The wave functional can also be expressed in momentum space, 

1 f d 3 k 



¥[ 



/Vexp 



(2tt) 3 



i(k)0x(k) + 03(k)0a(k) + c0!(k)0 2 (k) + c03(k)0i(k)] ^(77) 



(22) 



where 



* fc = U k cxp ^ -- [^(k)^(k) + 0j(k)0 2 (k) + c0!(k)0 2 (k) + C0^(k)0i(k)] A fc fa) 



(23) 



and 0j(— k) = 0*(k) for the real scalar field. Because there is no coupling between modes with different k, we will 
only consider a single wavelength and drop the index k for convenience from now on. 



III. THE DENSITY MATRIX AND THE COHERENCE LENGTH 



We now have the wave functional for all modes with single wavelength k. The next step is to calculate the reduced 
density matrix for 0i by tracing out 2 . 



p(0i,0i;f?) = / d<j) 2 d4>* 2 ¥£(0i, 2 , 77)^(01, 2 , 77) 

1 



Wk\ 2 / d(t> 2 d(/)* 2 exp 



101 + 2 2 + C010 2 + C0 2 0i)A* 



10* + 020 2 + C010 2 + C0 2 0i)A 



(24) 
(25) 
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This can be computed from the Gaussian integral: 



where 



p((f>i,<Pi;v) 



R=-^ L (\cbi\ 2 +\^n + 



4n 



A + A 



-\Af k \ 2 exp(R + iI) 



(26) 



8(A + A*) 



i(A+A*) 2 i(\<t>i\ 2 + \ti\ 2 + <f>:ti + <f>i4>*i) 



il = -{l-c 2 ) 



+(A* - A) 2 (|0i| 2 + \4>i\ 2 - - <j>i4>t)] 
^f^(l0i| 2 -|^i| 2 ) 



To determine the coherence length of the reduced density matrix, it is convenient to introduce the new variables 

1 



X 



+ <h) 



l 



A = - fa -fa) 



In terms of these variables, the reduced density matrix (26) becomes 

2 



p(<t>i,<f>i;v) = A + A * l-^fel ex P 



X " + ^+/3( X A*+X*A) 



(27) 
(28) 

(29) 
(30) 

(31) 



Because /3 



-(A* — A) is purely imaginary, the third term in the exponential just gives a complex phase. The 



first term gives the dispersion of the system, the dispersion coefficient a being 



(l-c 2 ){A + A*) 



(32) 



The second term describes how fast the density matrix decays when considering the off-diagonal terms. Hence, l c is 
called the coherence length and is given by 



i, 



\ 



(A + A*) 



\-6 



(i+A A *) 



(33) 



IV. DECOHERENCE IN THE USUAL INFLATION MODEL 

f 



For usual inflation, aft) = e Ht which is equivalent to air/) = Here, H is the Hubble constant. Eq. (|21| then 

tells us 



-ikrj 



2k 



1 



ki] 



- C 2 



f 



(34) 



Considering the wave functional (23), we have to require a positive real part of A for obvious reasons. Therefore, we 
choose ci — and 

k 



Mv) = 



1 



#V i 



kit 



Then, Eq. ( 33 ) gives us the coherence length 1 : 



l,= 



H(l + k 2 V 2 ) 1 / 2 



k 3/2 (l + 



1/2 



(35) 



(36) 



We recover the results in I21| after accounting for some typos in that paper. 
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TABLE I: Table (comparing power law inflation and ekpyrosis) 





power law inflation 


ekpyrotic phase 


range ot t 


u s I \ oo 


— CO \ /, \ u 


a(t) 




(-ty 


P 


p » i 


p < i 


range of r\ 


— oo < r\ < 


—00 < T) < 


a(j]) 


[(l-p)^/(l-f) 


[-(l-p)i7]' /(1 - rt 


a 


p i 


p i 


a 


(i— p) V 


(l-p) 17 


a 


p(2p-l) 1 


p(2p-l) 1 


a 


(1-p) 2 >7 2 


(1-P) 2 »7 2 



We see that if no interaction is present (c = 0), the coherence length approaches a constant value. Adding even a 
small interaction will reduce it to zero (See Fig. [I]). Besides, the coherence length starts to decrease exponentially 
when the wavelength crosses the Hubble radius, which justifies our heuristic derivation in cosmological perturbation 
theory. 




FIG. 1: The relation of coherence length and the conformal time for usual inflation. The horizontal axis is kr\ and the vertical 
axis is normalized coherence length. The upper (red) line corresponds to no interaction, and the lower (blue) line corresponds 
to c = 0.15. If there is an interaction, the coherence length starts decreasing and eventually becomes zero for the superhorizon 
modes. 



V. DECOHERENCE IN POWER LAW INFLATION AND EKPYROTIC PHASE 

The scale factor behaviors of power law inflation and ekpyrosis are very similar so we consider them at the same 
time. We list some properties of their scale factors in the Table |T| 



Because both of the power law inflation and ekpyrosis have the same — , they share the same solution of Uk- The 

a 



differential equation of (21 1 can be solved exactly by 

u k = v 73 ^ [ciH^i-krf) + c 2 H^(-kri) 

(1 2) 

where Ha ' are Hankel functions, and we have defined 



a 1 

a 4 



1 - 3p 



2(1 -P) 



(37) 



(38) 
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As before, we want Ak(rf) to have a positive real part, so we take c\ = 0, and Eq. (20) tells us 

l-3p 1 k H^i-kv) - H^i-kri]' 



2(1 -p)v 2 



Hg\-kr,) 



(39) 



Notice that they are the same for both power law inflation and ekpyrotic phase except p> 1 for the former and p< 1 
for the latter. We can then use Eq. (33) to calculate the coherence length for both cases. The numerical solutions are 
plotted in the Fig. [2] and Fig. [3] 



k 

2.0 r 




FIG. 2: The relation of coherence length and the conformal time for power law inflation. We choose p — 10 in this plot. The 
upper (red) line corresponds to no interaction, and the lower (blue) line corresponds to c = 0.15. 




k,, 



FIG. 3: The relation of coherence length and the conformal time for ekpyrosis with p = 0.1. The upper (red) line corresponds 
to no interaction, and the lower (blue) line corresponds to c = 0.15. It is clear that even the modes go outside the horizon, the 
coherence length continues growing and approaches to a nonzero constant in the end. 



In order to get the behavior of the coherence length l c when the modes are well outside the Hubble radius, we need 
the asymptotic form of the Hankcl function as x — > 0: 

F(a) fx\~ a r(a— 1) /a;\ 2 - Q l 



1 



— (-) 



a +2 



+ i 



_r(a + 1) \2J r(a 

where a > and T(a) is the Euler gamma function. After some manipulation of algebra, we have 

2-2cT 

2 l-2 a | 1 _ p |l-2a fe 2a ' 



Mv) 



2 l-2a| 1 _ p |l-2 afc 2a 



1 



r(«) 2 

- 



a — 1 

„2 



-kr) 



+ i 



T(a) 2 2aT(a) 4 



—kr] 



2<v 



if a > I 



if a < 1 



(40) 



(41) 
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as — ki] <C 1. 

For power law inflation, p> 1, we have a = | + ^zy 



e, < e <C 1. Therefore, 



lo 



1 + c 2 



r(<*) 4 



(a - 1) 2 tt 2 



(42) 



where 



^ = |2(l-p)| 



2+2e; -3-2e 



r(«) 2 



(43) 



From Eq. (42), it is obvious that if no interaction is present, the coherence length approaches a constant value Iq. 



However, even a small interaction will reduce the coherence length to zero just like what happened in the usual 
inflationary case. 



As for the ekpyrotic phase, p <C 1, and a — \ — 



| - e, < e < 1. Use Eq. Ell), it is not difficult to get 



l r ~ 



1 + c 2 



4a 2 r(a) 4 V 2 



-krj 



2-4e 



(44) 



This means the coherence length approaches a nonzero constant value no matter whether the interaction is present 
or not, in agreement with our numerical results in Fig. [3} 



VI. CONCLUSION 



We have studied a simple model with two free scalar fields interacting via a gradient coupling term in three different 
background spacetime: the usual inflation, the power law inflation, and the ekpyrosis. We also calculate the reduced 
density matrix and the corresponding coherence length by summing over one of the fields in all three cases. 

Our results are that if no interaction is present, the coherence length approaches a constant value. Adding even a 
small interation will reduce it to zero in either usual inflation or power law inflation case. Since this decoherence starts 
at Hubble crossing, the quantum fluctuations evaluated at krj = — 1 give the classical initial density perturbations 
which become the seeds of inhomogenities of our universe later on. However, this argument does not work for 
ekpyrosis whose coherence length never hits zero. This means the quantum coherence would not disappear even when 
the modes leave the horizon. Therefore, the heuristic argument that the quantum fluctuation can become classical 
for supcrhorizon modes is not valid for ekpyrotic phase. The implication of our result is that the power spectrum 
of CMB fluctuations is not directly related to the ekpyrotic phase. Even though at the end of ekpyrosis the scalar 
field has a scale-invariant power spectrum, it is hard to say anything about what we observe right now, since that 
depends on the "classical" initial density perturbations. This puts some doubts on the analyses of the cosmological 
perturbations in the cyclic/ekpyrotic universe. 

We derived our results using a very simple model. In principle, if we would like to claim the decoherence phenomenon 
cannot occur in ekpyrosis, we have to consider all kinds of interactions between systems and environment which is 
almost impossible to do. However, we believe the basic physics should emerge from simple models. We can easily 
generalize our analyses to a massive scalar field, and the results wouldn't change too much. We could also consider 
different kinds of interactions, but we will leave it to the future work. 

Finally, we model the environment with a scalar field, which is convincing but might be an oversimplified assumption. 
The environment can also be taken to consist of the short wavelength modes which are coupled to the long wavelength 
modes via non-linear couplings [22-28 . Hence, this might be another possible way to generate decoherence during 
ekpyrosis. 
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